Weakly disordered absorbing-state phase transitions 
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The effects of quenched disorder on nonequilibrium phase transitions in the directed percola- 
tion universality class are revisited. Using a strong-disorder energy-space renormalization group 
method, it is shown that for any amount of disorder the critical behavior is controlled by an infinite- 
randomness fixed point in the same universality class of the random transverse-field Ising models. 
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Using the formalism and the knowledge of equilibrium 
phase transitions, a natural trend with the aim of es- 
tablishing and classifying possible universality classes in 
nonequilibrium transitions arose [1, 2]. It was conjec- 
tured that the critical behavior of short-ranged interact- 
ing models with scalar order parameter and absence of 
conservation laws and extra symmetries are in the Di- 
rected Percolation universality class [3, 4, 5], which sep- 
arates an active fluctuating state from an inactive (ab- 
sorbing) nonfluctuating one [6]. Examples include tran- 
sitions in the contact process [7], catalytic reactions [8], 
depinning interface growth [9, 10], and marginal growth 
of turbulent domains in laminar flows [11]. 

Despite the theoretical understanding on the ubiqui- 
tous Directed Percolation universality class, its critical 
exponents have hardly been seen in real experiments 
[12] (see, however, Ref. [13]). It was then suspected 
that quenched disorder may be responsible. For spa- 
tial dimension d < 4, this is indeed the case as dic- 
tated by the Harris criterion [14, 15, 16] and confirmed 
by field-theoretical methods [17], which showed that the 
renormalization-group equations have only runaway so- 
lutions towards large disorder. In addition, disorder- 
dependent Griffiths-like phases [18, 19] nearby criticality 
have been observed [16, 20, 21, 22, 23, 24, 25]. 

This scenario thus points out an unconventional crit- 
ical behavior originating from the interplay between 
large spatial disorder fluctuations and strong correla- 
tions. Motivated by this reasoning, a strong-disorder 
renormalization-group (SDRG) method [26, 27] was ap- 
plied to the random contact process model [28]. For 
strong disorder, the critical behavior is governed by a 
universal infinite-randomness fixed point (IRFP) in the 
same universality class of the random transverse-field 
Ising model, whose dynamical scaling is known to be ac- 
tivated [29, 30, 31], i.e., length £ and time r are related 
through bar ~ with ip (dubbed tunneling exponent) 
being universal. For weak disorder, on the other hand, 
the critical point has finite disorder and usual power-law 
scaling r ~ £ 2 with nonuniversal dynamical exponent 
z proportional to the disorder strength and is formally 
infinite at the transition between the weak- and strong- 
disorder limits. These conclusions were also supported 
by density-matrix renormalization-group calculations in 
d = 1 [28]. Further Monte-Carlo calculations in d = 2 



confirmed the above scenario. However, the possibility 
that the weak-disorder regime was an artifact of finite- 
size effects was raised [32]. 

Facing the logarithmically slow dynamics, large-scale 
Monte-Carlo simulations in d = 1 for system sizes up 
to 10 7 sites and times up to 10 9 were performed [33]. 
No nonuniversal weak-disorder critical regime was found, 
shrinking considerably the parameter space in which it 
would exist and, together with the field-theoretical re- 
sults, strongly suggesting its nonexistence. It then raises 
the following puzzle. How can the SDRG suggest a fi- 
nite-disordered fixed point while Monte-Carlo simula- 
tions point to an m/imie-disordered one? Since the 
SDRG method is devised to include any minimal effects 
of disorder, it should be able to capture the physics of 
any IRFP as well as to point out its existence. 

This Brief Report is devoted to solve this question. In 
generalizing the SDRG method, we show that the crit- 
ical system is governed by a universal IRFP when any 
amount of disorder is present. Moreover, our motivation 
goes beyond the issue of settling the correct universality 
class of weakly disordered absorbing-state phase transi- 
tions. It deals with the delicate issue of implementing 
a SDRG in such a limit, which is an important tool to 
tackle many disordered systems. 

For definiteness, we now introduce the system, review 
the usual SDRG for random contact process [28], point 
out its failure, and modify it in order to overcome this 
problem. 

The contact process can be defined in a lattice in which 
each site i can have either a healed (<r, = 1) or an in- 
fected {(Ji = — 1) particle. A healed particle at site i 
can be contaminated by an infected one in a neighboring 
site j at rate Xij = Xji [34]. Also, an infected particle 
at site i can get spontaneously healed at rate /i,. The 
system has a stochastic dynamics governed by a mas- 
ter equation dtP({cr},t) = — HP({a},t) where the vec- 
tor P gives the probability of finding the configuration 
{a} = (<7i, o"2, . . . ) at time t and 



(i) 



is the generator of the Markov process [32, 35, 36]. Here, 

I , Q = 
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and restricts the sum to nearest neighbors only. 

The usefulness of this "quantum Hamiltonian formal- 
ism" comes from the fact that the steady state probability 
distribution P({er},i — ► oo) coincides with ground state 
of H and that the long-time relaxation properties are ob- 
tained from the low-lying spectrum of H. Although H 
is in general non-Hermitian, some standard methods can 
still be used. 

For the disordered case, Ajj and //j are random in- 
dependent variables distributed according to P (A) and 
R(fi), respectively. In this case, the low- lying spectrum 
of H can be reached by the following recipe (for simplic- 
ity, we focus on the d = 1 case) : (i) Search for the fastest 
("high-energy") scale in the system Q = max{Aj, /ij}, (ii) 
integrate out locally the corresponding mode, and (hi) 
renormalize the remaining degrees of freedom. Those 
steps are the basis of the SDRG method [37]. 

When (ii.a) O = A2, particles on sites 2 and 3 can be 
considered as one since they will be mostly in the same 
state, i.e., either both healed or both infected. Then, 
(iii.a) one treats H = A 2 (n 2 Q3 + Q 2 n3) exactly and 
Hi = fi2^h + as a perturbation. H n has two 

twofold multiplets. In the ground (excited) one, parti- 
cles 2 and 3 are in the same (opposite) state. Hi lifts the 
degeneracy of the ground multiplet, which corresponds 
to the effective healing rate jl of the particle cluster 2 
and 3. In second order of perturbation theory, one finds 
that Hi = jlM, with 
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When (ii.b) fl = fi2, the particle at site 2 can be 
considered as healed for all times. Hence, (iii.b) one 
treats Ho = fJ, 2 M 2 exactly and Hi = XiiniQi + Q1712) + 
Xiin-iQ-i + Qifrz) perturbatively. Ho has two fourfold 
multiplets. The ground (excited) one refers to particle 
2 healed (infected). Hi then lifts the degeneracy which 
corresponds to an effective infection rate A between par- 
ticles 1 and 3. In second order of perturbation theory 
Hi = X(niQ 3 + Qin 3 ), with 



A = nx\iX 2 /n2 , with K\ = 1 



(3) 



Once set the recursion relations (2) and (3), flow equa- 
tions for -P(A) and R{[±) can be constructed and the fixed- 
point distributions obtained [30, 32]. In principle, this 
give the long-time behavior of the system. The multi- 
plicative structure of Eqs. (2) and (3) is very important. 
Under these transformations, P(X) and -R(/-t) become in- 
definitely broad at criticality for any amount of disor- 
der as long as < K^x < 1 [30, 32, 38]. However, for 
> 1 the SDRG becomes inconsistent for weak dis- 
order because the renormalized couplings are typically 
bigger than the decimated ones. It is thus tempting to 
interpret this result as a runaway flow towards weak dis- 
order in odds with the field-theoretical results [17]. As 
we show below, this is not the case. The generation of a 
transition rate larger than the decimated ones is unphys- 
ical. The numerical prefactor > 1 is just an artifact 
of treating Hi until second order in perturbation theory. 
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Figure 1: (Color online) Schematic decimation procedure. 



According to Eq. (2), the splitting of the ground mul- 
tiplet of Ho due to Hi may overcome the distance (A2) 
between the two unperturbed multiplets for certain val- 
ues of fx 2 ,3 even though /x 2 ,3 < A2. Treating Ho + H\ 
exactly, however, this can never be the case. The ground 
state has energy and the excited ones are solutions of 
the polynomial 

X 3 - 2ix 2 + (^ + yU2M3) X ~ M2M3 (A 2 + £) = , (4) 

with £ = A2 + fi2 + M3- The renormalized healing rate jl 
is thus its minimal root. Although we could not solve jl 
analytically, its maximum value is shown to be /Umax = 
(2 — \/2)A 2 , which happens for ^3 = ^2 = A2. Numerical 
inspections of Eq. (4) show that jl < min{/j2, ^3, A2} in 
general. 

In addition, the operators connecting this particle clus- 
ter to the rest of the chain have also to be projected 
onto the same states. We find that 712,3 = &2,3n and 
Q2.3 = a2,3<2, where (1 + c 2 + 03)0-2,3 = 1 + c 2 ,3, with 
(A 2 + /J2,3-A) 2 C2,3 = A 2 /J3, 2 . (Note that 1/2 < a 2 ,3 < 1.) 
Therefore, the SDRG decimation procedure summarizes 
in replacing J2i=i,3 M n iQi+i+Qi n i+i)+J2i=2,3 ^i M i b y 
Xi(niQ+Qin)+jlQ+X3{nQ i +Qn i ), with Ai j3 = 0:2^1,3 
[see Fig. 1(a)]. The renormalization of Ai,3 is not consid- 
ered in the usual perturbative SDRG which is indeed a 
"weaker" effect since 012,3 € [1/2,1] and approaches 1 in 
the strong-disorder limit. 

Repeating the same procedure when decimating a heal- 
ing rate, Eq. (3) is then replaced by 



A = C-X 



(5) 



with 2( = Ai + A 2 + /i2 and x = \/C 2 A1A2, implying 
A < min{Ai, A 2 , /i 2 }- [Its maximal value A max = /j,2(3 — 
v5)/2 happens for Ai = A2 = /^ 2 -] Moreover, M\,3 = 

/3i,3-Mi,3, m,3 = n.1,3 and Qi, 3 = Qi,3, with i(3ix(C + 
\) = Ai(3/i 2 - A 2 ) + (A 2 + /J2X/U2 + A 2 + 2%) and /? 3 
is obtained by exchanging Ai ^ A2 in j3i. (Note that 
3/4 < Pi,3 < 1.) These results mean we have to replace 
J2i=o,3 X l (n l Q l+ i + Q l n i+ i) + Y Ji= i^^iM l by A (n Qi + 
Q ni) + piMx + X(hiQ 2 + Qin 2 ) + M3-M3 + A3(n 3 Q 4 + 
Q3^4), where fti^ = /?i,3/^i,3 [see Fig. 1(b)]. Note that 
"■1,3 (Qi,3) has no projection onto 713,1 (Q3,i). If this 
was not the case, the technical treatment of this SDRG 
would be more difficult because further-nearest-neighbor 
interactions would arise. Long-ranged interactions may 
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point out delocalized states. Their absence suggests that 
the SDRG here presented is amenable. 

Importantly, there are no level crossings in the entire 
region where the parameters of H± are less than or equal 
to the parameters of Hq, meaning the interpretation of 
the decimation steps still holds. Also important, the en- 
ergy difference between the second and first excited mul- 
tiplets (A21) of Hq + Hi only increases when increasing 
the perturbation and is always greater than the energy 
difference between the first excited and ground multiplets 
(Am). Precisely, A10 < (l - l/y/2) A 21 . 

Therefore, exactly projecting the entire system in the 
low-energy states of H + Hi makes the renormalization- 
group approach totally consistent. Whether or not these 
new recursion relations drive the system to the universal 
IRFP is not straightforwardly clear. This is the question 
we address in the next part of this paper. 

The fate of the critical point is obtained by solving 
the standard flow equations [30, 32] for P(A) and 
with the perturbed renormalized rates (2) and (3) re- 
placed by their exact counterparts (4) and (5) in addi- 
tion to the weaker renormalization of the neighboring 
transition rates (Ai,3 and fii^ in Fig. 1). Because of the 
complicated analytical structure of these quantities, a de- 
tailed analytical solution is hampered. We then rewrite 
j5 = K^/i2M3/A2 and A = n' x \i\2/ (12, where k'^ a are func- 
tions of the decimated transition rates. Moreover, we will 
neglect the renormalizations of A 13 and pi.3 [39]. Now, 
recall that (i) A and p, are always less than the decimated 
ones and that (ii) there is no correction to Hq in first or- 
der of perturbation theory. Point (i) permits us to set 
k'x = 1 in the weak-disorder limit. Hence, the system 
rapidly flows towards stronger disorder. As intermediate 
disorder is reached, the only way of stopping its further 
growth is making all decimations of type fl = const x /12 
[40], which corresponds to corrections in first order of 
perturbation theory. Point (ii) thus guarantees there is 
no hindrance on the flow towards even stronger disorder, 
in which limit n'^ x can be neglected [38] . We thus finally 
conclude that any amount of disorder drives the critical 
system towards the universal infinite-randomness fixed 
point. 

This conclusion was checked by numerical implemen- 
tation of the SDRG for weak- and moderate-disordered 
chains. (For consistency with the above proof, the weaker 
corrections to Ai j3 and p,i t 3 were neglected [41].) Follow- 
ing time r and length £ scales along the critical SDRG 
flow, the predicted [29] tunneling exponent V = 1/2 was 
confirmed (see Fig. 2). Here, r _1 = Q and is the den- 
sity of active particle clusters. The off-critical Griffiths 
phases surrounding the critical point in which r ~ £ 2 
with disorder-dependent dynamical exponent is also con- 
firmed in our numerics [42]. 

We now address the issue of weak disorder in higher 
dimensions. One key feature hinders a straightforward 
generalization of the RG steps here proposed: the recon- 
nection of the lattice. Because the coordination number 
increases in d > 1, one eventually needs to treat exactly 



big particle clusters. Leaving this task open, we cannot 
guarantee that all the RG steps will consistently lower 
the energy scale and drive the system to an IRFP. How- 
ever, and somewhat surprisingly, it was shown that the 
lattice reconnection does not hinder the flow towards in- 
finite randomness [31]. As the inconsistency of the sim- 
ple recursion relations is just an artifact of the perturba- 
tive treatment, it is then reasonable to conclude that the 
RG flows of the Directed Percolation and the Transverse- 
Field Ising universality classes are the same in the pres- 
ence of disorder in d — 2 and 3, as suspected in Ref. [32]. 
This would be in agreement with Monte-Carlo simula- 
tions in d = 2 [43], with the Harris criterion [14, 15, 16] 
and with the field-theoretical runaway flow solutions [17]. 

Recently, the experimental realization of the clean Di- 
rected Percolation universality class have been claimed 
[13, 44]. These experiments now raise another puzzle 
in the face of our results and many others [14, 15, 16, 
17, 32, 33, 43]. We would like to point out two crossovers 
which may give an explanation. One is the time crossover 
which was stressed in Ref. [33] (see, e.g., Fig. 7 therein). 
Because of the logarithmically slow dynamics, the "true" 
steady state takes place only after a long period of relax- 
ation. The other one is the clean-dirty crossover length. 
As in spin chains [45, 46], there is a crossover length 
below which disorder is irrelevant. The cleaner the sam- 
ple the longer the crossover length which reaches hun- 
dreds of sites even for spin chains with moderate disor- 
der. The time crossover is analogous to the tempera- 
ture crossover in spin chains. Only at very low temper- 
atures are the low-energy states important. The length 
crossover is equally analogous. Statistically rare fluctua- 
tions (the so-called large rare regions) only exist on large 
samples. Naturally, these crossovers are related through 
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Figure 2: (Color online) Time r and length £ scales along 
the SDRG flow at criticality. Transition rates are drawn from 
boxlike distributions as indicated. Chains have 2 x 10 6 sites 
and the data were averaged over 100 disorder realizations. 
Error bars are about the symbol size. 
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the dynamics. In Ref. [44], the system size is of order 
of hundreds of degrees of freedom. It is thus reasonable 
that the exponents measured are nonuniversal between 
the clean and the infinite-randomness fixed point. (This 
also may apply to other experiments [12].) The crossover 
length of the samples in Ref. [13] seems much bigger. 

In the face of the possibility of explaining many ex- 
periments, it is thus desirable to study the aforemen- 
tioned crossover of the exponents, which should be ac- 
complished without much effort by Monte-Carlo calcu- 
lations in d — 1, for instance. From the experimental 
side, it is desirable to pinpoint precisely the source of 
quenched disorder and to estimate its strength. Border 
effects may also diminish the effective size of the sam- 
ple. Finally, due to the slow relaxation processes, time 
measurements have to be carefully taken when locating 
the critical point. These studies should shed considerable 
light on this problem. 

In conclusion, we have modified the usual strong- 



disorder renormalization-group method in order to ex- 
actly recast the low-energy spectrum of the local fast- 
mode Hamiltonian. This allowed the method amenable 
to attack the problem in the weak-disorder limit in which 
the perturbative treatment yielded to runaway flow to- 
wards weak disorder. Applications to quantum spin 
chains as well as comparison with similar generalizations 
will be presented elsewhere. 

As discussed in Refs. [30, 38], this renormalization 
group method is not justified in the weak-disorder limit. 
We, however, leave open the possibility that, by exactly 
projecting the entire Hamiltonian onto the local low- 
energy spectrum, the method will correctly point out 
whether weak disorder is irrelevant. 

We are indebit to T. Vojta, E. Miranda, M.-Y. Lee, 
and K. A. Takeuchi for useful discussions. This work was 
supported by the NSF under Grants Nos. DMR-0339147 
and DMR-0506953, and by Research Corporation. 
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